Nonlinear spinor field in Bianchi type-I cosmology: accelerated regimes 
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A self-consistent system of interacting nonlinear spinor and scalar fields within the scope 
of a Bianchi type-I cosmological model filled with perfect fluid is considered. Exact self- 
\^ \ consistent solutions to the corresponding field equations are obtained. The role of spinor 

O I field in the evolution of the Universe is studied. It is shown that the spinor field gives rise to 

^ ' an accelerated mode of expansion of the Universe. At the early stage of evolution the spinor 



field nonlinearity generates the acceleration while at the later stage it is done by the nonzero 



^ , spinor mass. 

< 
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I. INTRODUCTION 



The accelerated mode of expansion of the present day Universe encourages many researchers 
O ■ to introduce different kind of sources that is able to explain this. Among them most popular is the 
dark energy given by a A term |[lll2.0]> quintessence |0>BS.0]> Chaplygin gas fs","^. Recently 
^1 cosmological models with spinor field have been extensively studied by a number of authors in a 
series of papers [10, The principal motive of the papers [ 

was to find out the regular solutions of the corresponding field equations. In some special cases, 
rS \ namely with a cosmological constant (A term) that plays the role of an additional gravitation field, 
■ we indeed find singularity-free solutions. It was also found that the introduction of nonlinear 
spinor field results in a rapid growth of the Universe. This allows us to consider the spinor field 
as a possible candidate to explain the accelerated mode of expansion. Note that similar attempt is 
made in a recent paper by Kremer et. al. fllQL In this paper we study the role of a spinor field 
in generating an accelerated mode of expansion of the Universe. Since similar systems, though 
from different aspects were thoroughly studied in [ 13, 14], to avoid lengthy calculations regarding 
spinor and scalar fields, we mainly confine ourselves to the study of master equation describing 
the evolution of BI Universe. We here give the solutions to the spinor and scalar field equations, 
details of these solutions can be found in [il3i.il4]. 



n. BASIC EQUATIONS : A BRIEF JOURNEY 

We consider a self consistent system of nonlinear spinor and scalar fields within the scope of a 
Bianchi type-I gravitational field filled with a perfect fluid. The spinor and the scalar field is given 
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by the Lagrangian 

^ 2 



mVirvA + F + ^(l+AiFi)(p,«<, (2.1) 



where X\ is the coupling constant and F and F\ are some arbitrary functions of invariants generated 
from the real bilinear forms of a spinor field. Here we assume F — F{I,J) and Fi = Fi{I,J) with 
I = S^, S= xffXj/, J = P^ andP = ify^Xj/. 
The gravitational field is chosen in the form 

ds^ = dt^ — a^dxl — a\dx\ — a^dx^, (2.2) 

where a,- are the functions of t only and the speed of light is taken to be unity. We also define 

T = aia2a2. (2.3) 

We consider the spinor and scalar field to be space independent. In that case for the spinor and 
the scalar fields and metric functions we find the following expressions |fl4ll. 

For F = F{I) we find S = Cq/t with Cq being an integration constant. The components of the 
spinor field in this case read 

¥i,2{t) = (Cl,2/v^)e-'^ V/3,4(0 = (C3,4/v^)e'^ (2.4) 

with the integration constants obeying Cq as Cq = C\-\-C\ — Cl —Cl- Here /3 = j{m — &)dt with 
& = dF/dS + {Xi(p^/2)dFi/dS. 

For F = F{J) in case of massless spinor field we find P — Dq/x. The corresponding compo- 
nents of the spinor field in this case read: with 



Wi,2 = {Di2e''' + iD3Ae-"' 

¥3.4 = {iD,y+D3,4e-''')/V^, 

with Do = 2 {Dj +Dl-Dl- Dj) . Here o = j'^dt with = dF jdP + ( Ai ^'^/2)dFi /dP. 
For the scalar field we find 

dt 



(2.5) 



where C and Ci are the integration constants. 

Solving the Einstein equation for the metric functions we find 

a,{t) =A,[x{t)Yl\x^[X,j [x{t')]-Ut'], (2.7) 

with the integration constants A/ and Xi obeying A1A2A3 = 1 and Xi +X2 +X3 = 0. Note that 
to evaluate the metric functions at any given time t we should first integrate / and only then 
substitute t by t. 

The theoretical arguments and recent experimental data which support the existence of 
an anisotropic phase that approaches an isotropic one, led us to consider the models of Universe 
with anisotropic background. On the other hand the isotropy of the present-day Universe lead 
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us to study how the initially anisotropic BI space-time can evolve into an isotropic Friedman- 
Robertson-Walker (FRW) one. Since for the FRW Universe ai{t) = a2{t) = a3{t), for the BI 
universe to evolve into a FRW one we should set Di = D2 = = I. Moreover, the isotropic 
nature of the present Universe leads to the fact that the three other constants X, should be close to 
zero as well, i.e., << 1, (/ = 1,2,3), so thatXif[T{t)]^^dt for f < 00 (for T{t) = t" with 
n> \ the integral tends to zero as ? — > 0° for any Xj). The rapid growth of the Universe due to the 
introduction of the nonlinear spinor field to the system results in the earlier isotropization. 

As is seen from (12.41 1. ( 12.51) . (12.61 ) and ( 12.71) . the spinor, scalar and metric functions are in some 
functional dependence of T. It should be noted that besides these, other physical quan tities such as 
spin-current, charge etc. and invariant of space-time are too expressed via T lf l3Lll4ll . It should be 
noted that at any space-time points where t = the spinor, scalar and gravitational fields become 
infinity, hence the space-time becomes singular at this point [ 14]. So it is very important to study 
the equation for T (which can be viewed as master equation) in details, exactly what we shall do in 
the section to follow. In doing so we analyze the role of spinor field in the character of evolution. 



III. EVOLUTION OF BI UNIVERSE AND ROLE OF SPINOR FIELD 

In this section we study the role of spinor field in the evolution of the Universe. But first of all 
let me qualitatively show the differences that occur at the later stage of expansion depending on 
how the sources of the gravitational field were introduced in the system. In doing so we write the 
Einstein equation in the following form: 

0.2 ^ ^ CI2 ci^ 
a2 as a2 

a^i a\ gt, a\ 
d\ ^0,2 d2 
ai CI2 cii CI2 
d\ d2 ^^2^3^ dj, d[ 

a\ CI2 ^2 ^3 CI3 ci\ 

Here A is the cosmological constant, is the energy-momentum tensor of the source field. The 
Eq. (13.11 ) is thoroughly studied in [13]. After a little manipulation from (13.11 ) one finds the equation 
for T which is indeed the acceleration equation and has the following general form: 

^- = \k(tI + T^')+3A, (3.2) 

Note also that here a positive A corresponds to the universal repulsive force which is often consid- 
ered as a form of dark energy, while a negative one gives an additional gravitational force. 
The Bianchi identity G^.^ = in our case gives 

t^--l{T^-Tl). (3.3) 





(3.1a) 


KTi+A, 


(3.1b) 




(3.1c) 


kT^+A. 


(3.1d) 



After a little manipulation from (13.21) and (13.31) one finds the following expression for Tq-. 

kT^ = 3H^~A-Coo/t\ (3.4) 
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where the definition of the generalized Hubble constant H as 




(3.5) 



Let us analyze the relation (13.41) in details. Consider the case when A = 0. At the moment when 
the expansion rate is zero (it might be at a time prior to the "Big Bang", or sometimes in the far 
future when the universe cease to expand we have H = 0.) the nonnegativity o f suggests that 
Coo < 0. Before considering the case for large T we should like to study the Eq. (13. J I) in detail. For 
the spinor and scalar fields chosen in this paper they are identically fulfilled. If this is not the case, 
an additional equation, know as equation of state, is applied to connect pressure (T^) with energy 
density (Tq). In the long run from (13.31) one finds something like (7o')*T = const., where b is some 
constant (in case of perfect fluid b = I + Q. Thus we see that the energy density of the source 
field introduced into the system as above decreases with the growth of T. Now if we consider the 
case when t is big enough for Tq to be neglected, from (13.41) we find 



From (13.61) and (13.71) it follows that for T to be infinitely large, A > 0. In case of A = we find that 
beginning from some value of t the rate of expansion of the Universe becomes trivial, that is the 
universe does not expand with time. Whereas, for A > the expansion process continues forever. 
As far as negative A is concerned, its presence imposes some restriction on T, namely, there exists 
some upper limit for T (note that T is essentially nonnegative, i.e. bound from below). Thus we 
see that a negative A, depending on the choice of parameters can give rise to an oscillatory mode 
of expansion [13J. Thus we can conclude the following : 

Let be the source of the Einstein field equation; Tq is the energy density and r/, Tj^, 

are the principal pressure and = T2 = T^. An ever-expanding BI Universe may be obtained if 
and only if the A term is positive (describes a repulsive force and can be viewed as a form of dark 
energy) and is introduced into the system as in (13.11 ) or if the source field introduced as a part of 
energy -momentum tensor behaves like a A term as X —>■ 0°. 

It should be noted that the sources of the gravitational field such as spinor, scalar and electro- 
magnetic fields, perfect or imperfect fluids, as well as dark energy such as quintessence, Chaply- 
gin gas are introduced into the system as parts of the total energy-momentum tensor T^. It is also 
known that the dark energy was introduced into the system to explain the late time acceleration of 
the Universe. To show that though the dark energy is introduced into the system as a part of total 
energy-momentum tensor, it still behaves like a A term as x ^ 0°, we write them explicitly. The 
quintessence and Chaplygin gas are given by the following equation of states: 



3//2-A^0. 



(3.6) 



On account of (13.51 ) from ( 13.61 ) one finds 




(3.7) 



Pq = W£q, W e [-1,0], 

Pc = -A/e,., A>0. 



(3.8a) 
(3.8b) 



Note that the energy densities of the quintessence and Chaplygin gas are related to T as 



l+w 



WG[-1,0], 



(3.9a) 




(3.9b) 
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From i3.% and (13 .81) follows that £c \fA and pc — > — Va as t ^ oo. In case of a quintessence, 
for w > —1, both energy density and pressure tend to zero as T tends to infinity. But for w = — 1 
(sometimes known as phantom matter) we have £^ Eoq and pq —£oq as T ^ <». It means a 
quintessence with w = — 1 and Chaplygin gas behave like a A term when r ^ oa and hence can 
give rise to an ever expanding Universe. 

Before solving the equation for T we have to write the components of the energy-momentum 
tensor of the source fields in details: 

= mS-F + ^{l+XiFi)(p^ + £pf, 

(3.10) 

Ti = Ti = Ti = ^S + ^P-F-^-{l+XiFi)(p^-ppf, 

where, ^ = 2SdF /dl + XiScp^dFi/dl and ^ = 2PdF /dJ + liPcp^dFi/dJ. In dTTOb £pf and ppf 
are the energy density and pressure of the perfect fluid, respectively and related by the equation of 
state 

Ppf = C£pf, Ce[0,l]. (3.11) 

Let us now study the equation for T in details and clarify the role of material field in the 
evolution of the Universe. For simplicity we consider the case when both F and Fi are the functions 
of / (5) only. We also set C = 1 and Cq = I. Thanks to the spinor field equations and those for the 
invariants of the bilinear spinor form, the energy-momentum conservation law for the spinor field 
satisfied identically lfl3l] . As a result the Eq. (13.31) now reads 1 13i] 

£ + -{£+p)=0. (3.12) 

In view of (13.1 II) from (13.121) for the energy density and pressure of the perfect fluid one finds 

_ gp _ Co^o 

Further we set £o = 1. Assume that F = XS'^ and Fi = S'' where A is the self-coupling constant. 
As it was shown in lUsll . the spinor field equation, more precisely the equations for bilinear spinor 
forms, in this case gives S = Cq/t. Then setting Cq = 1 for the energy density and the pressure 
from dTTOb we find 

m A t'~-^ 1 

(3.13) 

(g-l)A [{\-r)Xi + T'-]T'-^ g 

2(Ai + T'-)^ t1+^~^' 

Taking into account that Tq and are the functions of T only, the Eq. (13.21) can now be 
presented as 

T = ^{qi,T), (3.14) 

where we define 




^(^i,T) = (3/2)fc(m + A(^-2)Ti-^ + AirT'-V2(Ai + T'')V(l-C)/T^), (3-15) 
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where qi = {fc, m, A, Ai,^, r, ^} is the set of problem parameters. The En. (13.141) allows the fol- 
lowing first integral: 

i=^/2[E-^{q,,T)] (3.16) 

where we denote 



^(^l,T) = _-jK(^mT-A/T^-2_;Li/2(Ai + T'') + Tl-^ 



(3.17) 



From a mechanical point of view Eq. (13.141) can be interpreted as an equation of motion of a single 
particle with unit mass under the force ^{qi , t) . In (13.161) E is the integration constant which can 
be treated as energy level, and {q\, x) is the potential of the force ^{q\, x). We solve the Eq. 
( 13.141) numerically using Runge-Kutta method. The initial value of x is taken to be a reasonably 
small one, while the corresponding first derivative t is evaluated from (13.161) for a given E. 

Let us go back to the Eq. (13.141) . In view of (13.151 ) one sees, x (3/2) fcm > as T — > <», 
i.e., if X is considered to be the acceleration of the BI Universe, then the massive spinor field 
essentially can be viewed as a source for ever lasting acceleration. Note that it does not contradicts 
our previous statement about the role of energy-momentum tensor on ever expanding Universe, 
since the spinor field satisfies the Bianchi identity identically. 

Now a few words about considering f as acceleration. The Einstein equations for the FRW 
model read 



fa 

2-+ - 
a \a 



31-)' 



1 5 



From (13.181) one finds 



a 
a 



(3.18a) 
(3.18b) 

(3.19) 



The equation (13.191) is known as the acceleration equation. In analogy for the BI Universe from 
(13.11) we can write 



a] aj aj, 
a\ a2 CI3 



3Tl), 



(3.20) 



and declare it as acceleration equation. Though setting ai = = as we recover the original 
definition, hardly it will be helpful in our case. So in BI Universe we assume f be the acceleration 
and Eq. (13.21) be the acceleration equation. 

Let us now define the deceleration parameter. In FRW cosmology the deceleration parameter 
has the form 



aa 



1 + 



1 



(3.21) 



where //ftw = d/ais the Hubble parameter for FRW model. In analogy we can define a deceleration 
parameter as well. If we define the generalized deceleration parameter in the following way: 



d = - 



1 + 



H1+H2 + 
Hi + Hi + Hi 



(3.22) 



where Hi = di/at, then the standard deceleration parameter is recovered at a\ = a2 = aj,. But is 
this case the definition for acceleration adopted here is no longer valid. So we switch to the second 
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choice and following Belinchon and Harko et. al. 1 18 , [T^ define the generalized deceleration 
parameter as 

After a little manipulation in view of (13.11) and (13 .41) the deceleration parameter can be presented 
as 

d = -- ^ \ (3.24) 

Let us now go back to the equations (I3.14I) . (13.151) . (13.161) and (13.171) . As one sees, the positivity 
of the radical imposes some restriction on the value of T, namely in case of A > and q>2 the 
value of T cannot be too close to zero at any space-time point. It is clearly seen from the graphical 
view of the potential [cf. Fig. [l]|. Thus we can conclude that for some special choice of problem 
parameters the introduction of nonlinear spinor field given by a self-action provides singularity- 
free solutions. As it was shown in [13] the regular solution is obtained only at the expense of 
broken dominant-energy condition in the Hawking-Penrose theorem. 

If, in an eigentetrad ofT^y, e denotes the energy density and pi, p?, denote the three prin- 
cipal pressure, then the dominant energy condition can be written as ll20ll : 

e + I^Pa > 0; (3.25a) 

a 

e + Pa > 0, a = 1,2,3. (3.25b) 
The dominant energy condition for the BI metric can be written in the form: 



7-0 
''O 


> 




7-0 
^0 


> 


ii ai 


7-0 

''O 


> 


i2«2 


7-0 

''O 


> 


''3 <^3 



'3 ^^3' 



(3.26a) 
(3.26b) 
(3.26c) 
(3.26d) 

In Fig. |2lwe plot the potential for a negative A . As one sees, in the vicinity of T = there exists 
a bottomless potential hole. As one sees, if in case of a self-action the initial value of x is too close 
to zero and the constant E is less than '^max (the maximum value of the potential in presence of a 
self-action), the Universe will never come out of the hole. 

For numerical solutions we set K = 1, spinor mass m = 1, the power of nonlinearity we choose 
as ^ = 4, r = 4 and for perfect fluid we set ^ = 1/3 that corresponds to a radiation. We also set 
Coo = —0.001 and E = 10. The initial value of T is taken to be To = 0.4. The coupling constant is 
chosen to be Ai = 0.5, while the self coupling constant is taken to be either A = 0.5 or A = —0.5. 
Here, in the figures we use the following notations: 

1 corresponds to the case with self-action and interaction; 

2 corresponds to the case with self-action only; 

3 corresponds to the case with interaction only. 



As one sees from Fig. \\\ in presence of a self-action of the spinor field with a positive A, there 
occurs an infinitely high barrier as T ^ 0, it means that in the case considered here T cannot be 
trivial [if treated classically, the Universe cannot approach to a point unless it stays at an infinitely 
high energy level]. Thus, the nonlinearity of the spinor field provided by the self-action generates 
singularity-free evolution of the Universe. But, as was already mentioned, this regularity can be 
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FIG. 1: View of the potential '^(t) for A > 0. FIG. 2: View of the potential -^{z) for a nega- 

tive A. 



achieved only at the expense of dominant energy condition in Hawking-Penrose theorem. It is also 
clear that if the nonlinearity is induced by a scalar field, T may be trivial as well, thus givi ng ri se to 
space-time singularity. Note that cases in presence of a A term are thoroughly studied in lfl3l[l4ll . 
It was shown that introduction of a positive A just accelerates the speed of expansion, whereas, a 
negative A depending of the choice of E generates oscillatory or non-periodic mode of evolution. 
Note also that the regular solution obtained my means of a negative A in case of interaction does 
not result in broken dominant energy condition 1 14]. 

In Figs. |3land|4|we plot the corresponding energy density and pressure. As one sees, in case 
of a positive A the energy density is initially negative while the pressure is positive. In this case 
though the solution is singularity-free, the violation of dominant energy takes place. In case of a 
negative A the pressure is always negative. 

The purpose of plotting the energy density and pressure is to show that the energy density of 
the source field indeed decreases with the increase of the Universe. This also shows that there 
exists an interval where the energy density of the system with spinor field nonlinearity generated 
by the self-action is negative. This is in line with our earlier assumption. Moreover, we see the 
pressure of the source field becomes negative in course of evolution (In case of self-action with a 
positive X pressure is initially positive, but with the expansion of the Universe it becomes negative, 
whereas, in case of of a negative X as well as in case of interacting fields the pressure is always 
negative). Recall that the dark energy (e.g. quintessence, Chaplygin gas), modelled to explain the 
late time acceleration of the Universe, has the negative pressure. So we argue that the models with 
nonlinear spinor field and interacting spinor and scalar fields to some extent can be considered as 
an alternative to dark energy which is able to explain the late time acceleration of the Universe. 

In the Figs. |5]and|6lwe illustrate the acceleration of the Universe for positive and negative X, 
respectively. As one sees, in both cases we have decreasing acceleration that tends to (3/2) fcm as 

T — > °°. 

In Figs. |71and[8]we plot the deceleration parameter. 

The Figs. |5l|6l|7] and [H show the accelerated mode of expansion of the Universe. As one sees. 
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FIG. 3: Energy density and pressure correspond- 
ing to a positive A. 



FIG. 4: Energy density and pressure in case of a 
negative A. 
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FIG. 5: Acceleration of the Universe corre- 
sponding to a positive X. 



FIG. 6: Acceleration of the Universe in case of 
a negative A . 



the acceleration is decreasing with time. Depending of the choice of nonlinearity it undergoes an 
initial deceleration phase. It is also seen that the nonlinear term plays proactive role at the initial 
stage while at the later stage spinor mass is crucial for the accelerated mode of expansion. 
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FIG. 7: Deceleration parameter corresponding 
to a positive A. 



FIG. 8: Deceleration parameter in case of a neg- 
ative A . 



rV. CONCLUSION 

We considered a system of interaction nonlinear spinor and scalar fields within the scope of 
a BI cosmological model filled with perfect fluid. The spinor field nonlinearity gives rise to an 
effective negative pressure in the course of evolution. Comparison of the effective pressure of the 
nonlinear spinor field with that of a dark energy given by a quintessence or Chaplygin gas leads 
us to conclude that the spinor field can be seen as an alternative to the dark energy able to explain 
the acceleration of the Universe. It was shown that the nonlinear spinor term is proactive at the 
early stage of the evolution and essentially accelerates the process of evolution, while at the later 
stage of evolution the spinor mass holds the key. Given the fact that neutrino is described by 
the spinor field equation and it too possesses mass (though too small but nonzero), the presence 
of huge number of neutrino in the Universe can be seen as one of the possible factor of the late 
time acceleration of the Universe. It was also shown that for some specific choice of parameters 
it is possible to construct singularity-free model of the Universe, but this regularity results in the 
broken dominant energy condition of the Hawking-Penrose theorem. 



[1] Padmanabhan, T., Phys. Rep. 380 (2003) 235. 

[2] Sahni, V., Dark Matter and Dark Energy astro-ph/0403324. 

[3] Saha, Bijan, Anisotropic cosmological models with a perfect fluid and a A term (accepted for publica- 
tion in Astrophys. Space Sci.) gr-qc/0411080 

[4] Cladwell, R.R., Dave, R., and Steinhardt, RJ., Phys. Rev. Lett. 80 (1998) 1582. 

[5] Sahni, V. and Starobinsky, A.A., Int. J. Mod. Phys. D 9 (2000) 373. 

[6] Zlatev, I., Wang, L., and Steinhardt, PL, Phys. Rev. Lett. 82 (1999) 896. 

[7] Saha, Bijan, Anisotropic cosmological models with perfect fluid and dark energy reexamined (accepted 



Nonlinear spinor field in Bianchi type-I cosmology: accelerated regimes 



February 7, 20081 1 



for publication in Int. J. Theor. Phys.) gr-qc/0501067 

[8] Kamenshchik, A.Yu., Moschella, U., and Pasquier, V., Phys. Lett. B. 511 (2001) 265. 

[9] Saha, Bijan, Chinese J. Phys. 43 (2005), 1035-1043 

[10] Saha, B. and Shikin, G.N., J. Math. Phys. 38 (1997) 5305] 

[1 1] Saha, B. and Shikin, G.N., Gen. Relativ. Gravit. 29 (1997') IT099l 

[12] Saha, Bijan, Mod. Phys. Lett. A 16 (2001), 1287j 

[13] Saha, Bijan, Phys. Rev. D 64 (2001) 123501 

[14] Saha, Bijan and Boyadjiev, T., Phys. Rev. D 69 (2004) 124010*" 

[15] Armendariz-Picon, C. and Greene, PB., Gen. Relativ. Gravit. 35 (2003) 1637. 

[16] Ribas, M.O., Devecchi, F.P, and Kremer, G.M., Phys. Rev. D 72 (2005) 123502. 

[17] C.W. Misner, Astrophys. J. 151, 431 (1968). 

[18] J.A. Belinchon, Perfect fluid LRS Bianchi I with time varying constants, gr-qc/041 1005. 

[19] T. Harko and M.K. Mak, Int. J. Mod. Phys. D 11, 1171 (2002). 

[20] S.W. Hawking, and R. Penrose, Proc. R. Soc. London A 314, 528 (1970). 



